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ABSTRACT 

We use two independent methods to forecast the dark energy measurements achievable 
by combining future galaxy redshift surveys based on the radio HI emission line with 
Cosmic Microwave Background (CMB) data from the Planck satellite. In the first 
method we focus on the 'standard ruler' provided by the baryon acoustic oscillation 
(BAO) length scale. In the second method we utilize additional information encoded 
in the galaxy power spectrum including galaxy bias from velocity-space distortions 
and the growth of cosmic structure. We find that a radio synthesis array with about 
10 per cent of the collecting area of the Square Kilometre Array (SKA), equipped with 
a wide (10 — 100 deg 2 ) field-of-view, would have the capacity to perform a 20,000 deg 2 
redshift survey to a maximum redshift z max ~ 0.8 and thereby produce dark energy 
measurements that are competitive with surveys likely to be undertaken by optical 
telescopes around 2015. There would then be powerful arguments for adding collecting 
area to such a 'Phase-1' SKA because of the square-law scaling of survey speed with 
telescope sensitivity for HI surveys, compared to the linear scaling for optical redshift 
surveys. The full SKA telescope should, by performing a 20,000 deg 2 HI redshift survey 
to z max ~ 2 around 2020, yield an accurate measurement of cosmological parameters 
independent of CMB datasets. Combining CMB (Planck) and galaxy power spectrum 
(SKA) measurements will drive errors in the dark energy equation-of-state parameter 
w well below the 1 per cent level. The major systematic uncertainty in these forecasts 
is the lack of direct information about the mass function of high-rcdshift Hl-cmitting 
galaxies. 'Stacking experiments' with SKA pathfinders will play an important role in 
resolving this uncertainty. 

Key words: cosmology : radio surveys - galaxy power spectrum - baryonic oscillations 
- dark energy 



1 INTRODUCTION 



Over the past decade several different cosmological probes 
have collectively demonstrated that the Universe contains a 
large fraction of its energy density in clustered dark matter 
and unclustered dark energy. The analysis of the angular 
distribution of the Cosmic Microwave Background (CMB) 
anisotropics requires that the total energy density is close 
to that predicted in a spatially-flat Univers e. This data, 
given a strong prior on the Hubble consta nt dSpergel et al 



2007 ) imposed by the Hubble Key Project (|Freedman et al 



200 If ) , yields evidence for unclustered dark energy. In a Uni- 
verse that is close to spatially flat, observations of distant 
Type-la supernovae provide independent evidence for a dark 
energy component with a negative pressure driving an accel- 



erated expansion (|Riess et all 2006). Furthermore, the anal- 
ysis of the shape of the galaxy power spectr um at low red- 
shift requires a matter density fi m ~ 0.25 fe.g. lPercival et all 
120071 ). and the corresponding measurement of baryon acous- 
tic oscillations (BAO s) at z = 0.35 agrees with this picture 
jPercival et all 120071 ) . Other probes confirm this paradigm: 
a positive cross-correlation of the distribution of the photons 
of the CMB with the distribution of galaxies at low redshift 
indicates, via the Integrated Sachs- Wolfe (ISW) effect, that 
the growth of large-scale over-densities is consistent with 
the presence of dark energy with a negative pressure. The 
ISW effect has been detected at the > 2 . 5-cr level by several 
groups (e.g iBoughn fc Crittenden] 120051; ICabre et all [20061 ; 
iRassat et al.ll2006l ; iGiannantonio et alj|2006l ). 

No model has yet been successful in fitting all of the 
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above cosmological data without using a dark energy com- 
ponent. This does not rule out a modi fied gravity model (e.g. 
ISkordis et al]|2006l ; Im ena et ai1l2006l ) or more complex cos- 
mological assumptions s uch as back-reaction from inhomo- 
geneities (|Rasanenll2004l ) . although dark energy from back- 
react ion is disfavoured by other authors (jlshibashi fc Waldl 
120061 ). However, the addition of a dark energy component is 
the model that is the best-supported phenomenologically to 
date, even though we do not have a clear understanding of 
its origin or its nature. 

Although there is a wealth of information indicating 
the existence of a dark energy component in our Universe, 
the measurements of its properties, and specifically its equa- 
tion of state w — p/(pc 2 ) (where p and pc 2 are the pressure 
and energy density of dark energy respectively), are cur- 
rently poor. The simplest explanation for dark energy is Ein- 
stein's cosmological constant, parameterized by w = —1, al- 
though the observed magnitude of this vacuum energy poses 
a formidable difficulty for theory. Cur rent measurements of 
the parameter w are at the 10% level (|Spergel et al.ll2007l ). 

The discovery that w 7^ — 1 would yield vital informa- 
tion on the nature of dark energy. If w is a function of time 
(i.e. redshift) and space, it could potentially approximate to 
a constant close to —1 at late ti mes, depending on its na - 
ture and evolution with redshift (|Simpson fc Bridle! l200rj ). 
The dark energy would therefore behave as a cosmological 
constant at the present day. We need to measure both this 
equation-of-state parameter w and its evolution with time 
in order to prove or disprove whether dark energy is due to 
the presence of a cosmological constant. 

There are several types of observation which could 
in principle achieve a measurement of w with per-cent 
level accuracy, including weak gravitational lensing, super- 
novae, and cluster number counts. Future radio (HI emis- 
sion) redshift surveys might make significant contributions 
to this quest by measuring the large-scale distr ibution, and 
henc e clustering powe r spect r um, of galaxies ([Blake et al.l 



20041; iRawlings et all 12004 IXbdalla fc Rawlingsl I200I 



Zhang fc Pen! 120061 ; IChang et al.l 120081 ). The characteris 



tic length scale of 'wiggles' in the matter power spectrum 
due to BAOs provides a set of 'standard rulers' at any 
epoch (i.e. redshift) at which they can be detected. BAO 
tests are especially promising for measuring dark energy as 
they ar e believed to be relatively immune to systematic 
effects ( Blake fc Glazebro ok 2003|; ISeo fc Eisensteinl 120031 ; 
iHu fc Haimanll2003l ) . BAO measurements have already been 
successfully extr acted from the Sloan Digital Sky Survey 
(SDSS) dataset jEisenstein et al.ll2005l ; IPercival et al.ll20071 ) 
and, although the relatively small volume of the Universe 
probed by SDSS only permits a measurement of fi m rather 
than an accurate determination of w, the prospects look ex- 
tremely promising. 

Developing a strategy for performing redshift surveys in 
radio wavebands is very different from planning counterparts 
in optical light. We discuss in this paper the advantages and 
disadvantages of a radio survey compared to an optical sur- 
vey. We also make a prediction of how much information 
can be obtained by measuring the full galaxy power spec- 
trum, and compare this to an analysis retaining only the 
clean signature of the BAOs. 

In Sec[2]we discuss some of the differences between op- 
tical and radio surveys of large-scale structure. In Sec|3]we 



de tail how we simulate HI sur veys, also referring the reader 
to lAbdalla fc Rawlingsl (120051 ). In Secg]we present forecasts 
for measuring dark energy from several sizes of radio surveys 
that should, alongside CMB data from the Planck satel- 
be attainable over the coming decade. In Sec[5]we use 
Monte Carlo Markov Chain (MCMC) techniques to esti- 
mate the quality of future galaxy power spectra; see also 
lAbdalla fc Rawlingsl (|2007l ). In Sec|6] we comment on the 
mitigation of potential systematic errors in the analysis. We 
conclude our findings in Sec[7] 

Unless otherwise stated we adopt a ACDM model 
with an extra parameter w, a constant equation-of-state of 
dark energy. Ba sed on recent measu r ements of cosmologi - 
cal parameters (|Komatsu et al.l 120081 ; iDunklev et al.l 120081 ) 
we adopt the following choice of parameters for our fidu- 
cial cosmological model: {r2b,O c ,w,/i,n a ,crg,r} = {0.04,0.26,- 
1.0,0.72,1.0,0.9,0.09}, where Qt and Jl c are the present-day 
density of baryons and cold dark matter expressed as a frac- 
tion of the critical density, h = Ho/(100 km s" 1 Mpc" 1 ) 
represents the Hubble Constant, n s is the primordial spec- 
tral slope of the power spectrum that has an overall present- 
day normalization erg and r is the optical depth to the last- 
scattering surface at recombination. We consider only mod- 
els that are spatially flat (i.e. Oa = 0.7 in our fiducial cos- 
mology) and we ignore t he effects of neutrinos and other 
forms of hot dark matter (|Abdalla fc Rawlingsll2007l ). 



2 A CRITICAL COMPARISON OF OPTICAL 
AND RADIO REDSHIFT SURVEYS 

Undertaking 'all-sky' galaxy redshift surveys in optical and 
radio wavebands gives rise to very different issues. Although 
astronomers will be striving to complete both types of sur- 
vey over the coming decade, there is a natural competitive 
tension between the techniques that is discussed in this Sec- 
tion. 

The optical spectra of galaxies are much wealthier in 
information than their radio spectra. Optical galaxy spectra 
exhibit absorption features, and many have strong emission 
lines such as Lya, [Oil] and Ha, making redshift determi- 
nation straightforward with sufficient signal-to-noise. The 
information is sufficiently rich that one can pre-select a par- 
ticular redshift range or type of galaxy from the colours 
derived from a photometric imaging survey. The rest-frame 
equivalent widths of some emission lines can be very large for 
subsets of galaxies, e.g. the star-forming galaxies, so blind 
surveys are also plausible via 'slit-less' techniques. Care must 
be taken in both targeted and blind surveys to account for 
interlopers from emission-line mis-identifications. 

In the radio part of the spectrum the information in a 
galaxy spectrum is much more limited. The radio continuum 
is relatively featureless allowing for only the very crudest 
photometric redshift techniques, e.g. by using correlation s 
between redshift and radio spectral index (|Cruz et al.ll2007l ) . 
There is only one strong line in the radio spectrum - the 
'21-cm line' arising from hyperfine splitting in the hydrogen 
atom (HI) - and it is trivial for radio telescopes to obtain 
the spectral dimension as part of an imaging survey. Recent 
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technological breakthroughs allow this to be achieved over 
large frequency bandwidths. In principle, clean and blind 
redshift surveys are therefore straightforward, with the lu- 
minosity of the HI line providing a direct measure of the HI 
mass in an object. HI line surveys are hence similar to blind 
optical emission-line searches performed with Integral Field 
Units (IFUs) or slit-less spectroscopy. 

Unfortunately the HI line is extraordinarily weak. Thus, 
at the time of writing, all of the most important surveys 
of large-scale structure have been undertaken with optical 
telescopes. The state-of-the-art method of performing these 
surveys is to use wide-field (~ deg 2 ) multi-object spectro- 
graph that, because of the multiplex advantage of obtaining 
spectra of 10 2 ~ 3 objects per pointing, have been able to 
build up catalogues containing 10 5-6 redshifts using 10 3 ~ 4 
separate telescope pointings. With spectroscopic exposures 
of the order of hours these programmes take 10 2-3 clear 
dark observing nights and therefore many years to complete 
with a single telescope. Surveys such as the SDSS have now 
mapped out significant fractions (~ 1/4) of the sky, but only 
at low redshifts (z ^$ 0.5) where moderate-sized (4-m class) 
telescopes can measure absorption-line galaxy redshifts with 
relative ease. 

There are good prospects of pushing such studies to 
moderately-higher redshifts (z < 1) with 4-m-class ground- 
based telescopes. One technique is to target highly-biased 
tracers such as Luminous Red Galaxies, which require a 
lower target density to suppress shot noise in the power 
spectrum. In this manner the Baryon Oscillation Spectro- 
scopic Survey (BOSS)Q, an extension of the SDSS, will reach 
z ~ 0.7 over 25% of the sky by 2014. Another technique is 
to focus on objects that have strong emission lines such as 
[Oil], so that redshifts can be measured to a greater depth 
without requiring a detection of the galaxy continuum light . 
In this manner the WiggleZ survey ( Gl azebrook et alj |2007) 
at the Anglo- Australian Telescope (AAT) is probing galax- 
ies up to z ~ 1 over 2.5% of the sky by 2010. Both of these 
approaches require a method of isolating specific types of 
high- redshift galaxies, such as by using colour cuts and/or 
photometric redshifts. 

Optical surveys can be performed to still higher red- 
shifts using the greater light grasp of 8-m class ground- 
based optical telescopes and optimized observational strate- 
gies. The recent deployment on the 8-m Subar u Telescope of 
the F MOS near-IR multi-object spectrograph (|Dalton et al.l 
l2006t ) will enable large-scale projects at near-IR wave- 
lengths, allowing redshift surveys based on absorption-line 
(i.e. 4000-A break) systems at z <^ 1 and emission-line (i.e. 
Ha) redshifts up to z ~ 2. There are also plans for wide- 
field red-sensitive optical multi-object spectrographs such as 
the proposed Wide Field Multi-Object Spectrograph (WF- 
MOS); this instrument would possess > 1000 fibres and have 
sufficient sensitivity to push wide-field optical redshift sur- 
veys to z ~ 4. 

In principle, instruments like FMOS and WFMOS could 
survey the entire sky for a selection of galaxies. However, 
such surveys would be very slow for three main reasons: 
(i) pre-selection of targets is hard, as existing imaging sur- 
veys are either small or too shallow to select targets; (ii) the 



field-of-view (FOV) of these spectrographs is of order one 
square degree and exposures could take most of a clear dark 
night at z ^ 1.5 (with the exception of line emitting galax- 
ies); and (iii) surveys must share time with other astronom- 
ical programmes or instruments and are relatively limited 
in duration. In practice these future ground-based optical 
redshift surveys will probably be limited to < 2000 deg 2 of 
sky. Similar limitations on sky area apply to the alternative 
ground-based optical technique of obtaining 'blind' redshifts 
of high-redshift Lya emitters using multiple IFUs deployed 
in the focal plane, such as proposed by the HETDEX exper- 
iment (|Hill et aljboOSft . 

The most plausible approach for performing an all-sky 
high-redshift spectroscopic galaxy survey in optical or infra- 
red wavebands is to use a space-based platform, as sug- 
gested by several candidate proposals for the Joint Dark 
Energy Mission (JDEMJ0 funded in the United States or 
the 'Euclid' mission fRefregier fc the DUNE collaboration! 



120081 ; ICimatti fc the SPACE collaboration! 120081 ) funded m 
Europe. The extremely low sky background from space en- 
ables much faster surveys, and galaxy pre-selection can be 
avoided using slit-less techniques detecting high-redshift Hq- 
emitting galaxies up to z ~ 2. JDEM could be launched by 
2015 if funding continues. 

Future radio telescopes provide an alternative route to 
constructing all-sky galaxy redshift surveys up to z ~ 2. 
Radio telescopes benefit from a rapid rate of increase of sur- 
vey speed with collecting area. The characteristic observing 
time t map needed to detect all the HI line-emitting galaxies 
above a fixed line-flux limit over a fixed sky area scales with 
sensitivity and FOV in accordance with 

1 

Cmai 



oc fFOV 



(1) 



where / is the fractional sensitivity of the telescope, which 
we define relative to a full SKA (/ oc A e ff/T ays , where A e ff is 
the effective collecting area and T sys is the system temper- 
ature). The quadratic nature of Eq[T] implies a fundamen- 
tally different scaling to the corresponding etendue figure 
used to characterize the survey speed of optical telescopes 
(1/tmap oc A c hFOV in the case of optical redshift surveys). 
Future radio interferometers will possess collecting areas 
which are orders of magnitude higher than those in oper- 
ation today, and the technology exists to produce a field-of- 
view FOV ~ 100 deg 2 , significantly larger than achievable 
by optical spectrographs, and limited only by computational 
power. The large frequency bandwidth of the radio telescope 
will simultaneously map out a significant redshift range; the 
maximum redshift attained is determined by the sensitiv- 
ity of the radio telescope and the degree of concentration of 
interferometric baselines in the telescope core. 

In summary, although radio wavebands offer very lim- 
ited possibilities for filtering techniques to pinpoint distant 
galaxies similar to those commonly used in optical wave- 
bands, the sheer survey volume that future radio interfer- 
ometers can map suggests that radio telescopes have the 
potential to become very competitive facilities for produc- 
ing cosmologically-interesting redshift surveys. 

A significant uncertainty in modelling future radio sur- 
veys is that radio techniques have so far only yielded direct 



2 http://cosmology.lbl.gov/BOSS/ 3 http://jdcm.gsfc.nasa.gov/ 
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detections of HI em i ssion to z ~ . 2 (|Zwaan et al.l l200ll ; 
IVerheiien et alj|2007t ICatinella et alj|2008h . Various simple 
models extrapolating the p roperties of HI galaxies to h igher 
redshifts are considered by lAbdalla fc Rawlingsl (|2005l ); we 
emphasize here the critical role which SKA pathfinder tele- 
scopes will play in reducing these uncertainties to the level 
needed to opti mize the design of the final SKA. Pathfinders 
such as FAST (|Duffv et al.ll2008h . ASKAlfl and MeerKAffl 
will achieve significant expansions in frequency range and 
field-of-view over existing systems with similar sensitivity 
(i.e. A e ff/T ay8 ). One can then envisage a set of 'stacking 
experiments' that can yield statistical measures of HI at 
much higher redshifts than the z ~ 0.2 objects currently 
under study. The ~ 30 deg 2 FOV expected to be delivered 
by ASKAP will allow deep-field observations covering (in 
~ 3 pointings) ~ I0 4 'WiggleZ' star-forming galaxies up to 
z ~ 1. Stacking the values of radio flux density in channels 
corresponding to HI at the optical redshifts of these objects 
provides, in the absence of any systematic noise floor, an 
V 10 4 ~ 100-fold increase in effective sensitivity that more 
than compensates for the factor ~ 10 diminution of line flux 
expected between similar Hi-emitting galaxies at z ~ 0.2 
and z ~ 1. The sensitivity to be delivered by MeerKAT 
over smaller (~ 1 deg 2 ) sky patches will probe the general 
population of z ~ 1 galaxies: F MOS Guaranteed Time Ob- 
servations l|Dalton et al.1120061 ) expect to measure, via stel- 
lar absorption features as well as Ha emission, redshifts for 
~ 10 4 galaxies per square degree. These can also be stacked 
to achieve high-redshift statistical detections of HI, provided 
that the noise in the radio datasets is suitably well-behaved. 



3 SURVEY MODELLING 

We review here our method of simulating the HI surveys 
assumed in Sec|4] and S ecl5.2l In large part we h ave fol- 
lowed the prescription o f lAbdalla fc Rawlingsl d2005h which 
we now clarify in more detail (see also lZhang fc Per] l|200rJO . 



3.1 HI mass limit 

The r.m.s. sensitivity for a dual-polarization radio receiver 
at system temperature T ByB for an integration of duration t 
on a telescope of effective collecting area A e g is given by the 
usual radio telescope equatiorQ 



Srms — 



V2Ck B T s , 



A eS VAut 



1.6 /uJy 
~1~ 



Au/MRz V i/hr 



(2) 

where C i s a constant o f ord er unity (see Chap- 
ter 1 of iThompson et al.1 |200ll) . k B is Boltzmann's 
constant, Av is the frequency bandwidth, and / = 
[Aff/(10 6 m 2 )][r sys /50 K]" 1 . The value of / hence parame- 
terizes the sensitivity of the radio telescope compared to the 
SKA. 

In practice we used for our calculations the frequency- 
dependent telescopes sensitivity listed in Table 1.2 of 

4 http://askap.pbwiki.com/ 

5 http://www.ska.ac.za/ 



6 We note that the numerator of Eq.2 of lAbd alla & Ra wlingsl 
1120051 1 should read 100 mjy not 100 njy. 



iTavlor fc Braunl (|l999l ) which we express in the following 
fitting formula: 



(1.790365 v 2 - 4.53125 v + 4.46666) ^Jy 



30km/s /4hr 



AV V t 



(3) 



where v is the frequency of the redshifted HI line. The first 
term in Eq(3] is equivalent to a dual-polarization r.m.s. sen- 
sitivity in a velocity channel of width AV = 30 km s _1 in 
a t = 4 hr exposure. The underlying assumptions are that 
the scaling of S m a with telescope sensitivity / is indepen- 
dent of frequency, and C = V2, which makes the optimistic 
assumption that all statistically-independent signal samples 
can be successfully averaged over a given exposure time. The 
assumption that / is independent of frequency implies that 
the increasing contribution of sky temperature to T sys at 
lower frequencies is compensated for by an increase in A e g 
with decreasing frequency. Given that the brightness tem- 
perature of the Galactic foreground emission scales roughly 
as T Ga i ~ 20(^/408 MHz)- 2 7 K, in practice this will be- 
come increasingly difficult to realize below ~ 300 MHz (HI 
at z ~ 3.7) once the rapidly rising sky temperature begins 
to dominate the system temperature. 

The flux detection limit Sv im for galaxies is defined by 
the threshold parameter 



Tla — Sy lrn f Sri 



(4) 



We allow for a certain fraction of flux to be 'resolved 
out' using a model for the source sizes and the interfero- 
metric baseline distribution. The proper physical size of the 
sources is modelled as R(z) = Ro/(l + z) 2 following the dis- 
cussion in AR05, where Ro = 15 kpc. The fraction of flux 
detected from a source at redshift z (see also Fig. 2 of AR05) 
is then described as an integral over the distribution f(B) 
of baselines B: 



Frac detected = 



f(B) exp 



B 1 



2B {zf 



dB (5) 



where Bo = [Ax(z)]/[7ri?(,z)] and A = 0.21 m, x(z) is the 
co- moving radial distance to redshift 2, and the baseline 
distribution f(B) is given by 



/(*) = £• 



2 1 

~ R" eXP 

7T Bi 



_B 2 

'2B 2 



(6) 



where we use a sum of four different baseline components 
i = 1 to i = 4 where B 1 = 0.5 km, B2 = 3 km, B3 = 150 km 
and B4 = 2600 km (this latter corresponding to 0.02 arc- 
sec resolution for wavelength A — 0.21 m). The coefficients 
cii satisfy X^!=i a i = f (such that J f(B) dB = 1) and were 
taken as 01 = 0.2242, a 2 = 0.2746, a 3 = 0.2506, a 4 = 0.2506. 
This is a composite array and each en component of the ar- 
ray is distributed in a 'scale-free' configuration, as used in 
AR05. We note that the costs of data transport increase dra- 
matically for the longest baselines of a radio interferometer 
and thus the final realization of the SKA is likely to be more 
centrally concentrated than we have assumed, meaning that 
we are being conservative in our assumed ability of the SKA 
to detect HI in galaxies. 

The width of the HI emission line determines the rel- 
evant frequency bandwidth Av in Eq[2] and is related to a 



© 0000 RAS, MNRAS 000, 000-000 



Dark Energy with future HI surveys 5 



line-of-sight velocity width V(z) at redshift z: 

\ V(z) 



l + z 



(T) 



where vq = 1.4204 GHz is the rest-frame frequency of the HI 
line, c is the speed of light and V(z) is measured in the rest- 
frame of the galaxy. We follow AR05 in assuming a scaling 
with redshift for the typical rest-frame velocity width of the 
line V(z) = Vo/y/1 + z, where Vo = 300 km s _ , and ignore 
the effects of inclination. 

The luminosity of HI photons from a gas cloud of mass 
Mm can be derived using atomic physics (see, e.g., Sec. 2. 2 
of AR05). Assuming a flat profile for the emission line with 
frequency, we obtain the following equation linking Mm with 
the observed flux density S u : 

16ir mu 



Mb 



or in usual units: 



3 Ai2hc 



x{zY S v V(z) (1 + z) 



(8) 



M m = 0.235M Q (rr/Mpc) 2 (S„//xJy) (U/kms -1 ) (l+z) (9) 

where vhh is the atomic mass of hydrogen, A\2 is the Ein- 
stein coefficient for the transition, h is Planck's constant, 
and x(z) is the co-moving distance to redshift z. 




Figure 1. The HI mass function plotted at redshift 2 = 2 for a 
no-evolution model (black solid line), evolutionary models A (red 
dashed line), B (green dotted line) and C (blue dot-dashed line). 
This is a new version of Fig. 4 from AR05 which was incorrect due 
to missing factors of h in the plotting routine. The evolutionary 
models are described in Sec. 13.21 



3.2 HI mass function modelling 

We assume that the HI mass function at all redshifts is de- 
scribed by a Schechter function for the number density of 
galaxies dn in a mass range dM: 



An 
AM 



AM = 



M 
M* 



exp 



M 
M* 



AM 
~~M* 



(10) 



parameterized by a low-mass slope a, a characteristic mass 
M* and a normalization <f>* . This function has only been 
measured in the local Universe. The HIPASS survey re- 
ported the results: a = -1.3, M* = 3 .47 fo~ 2 10 9 Mq, 
</>* = 0.0204 h 3 Mpc -3 |Zwaan et all l2003l fl. Throughout 
this study we assume a — —1.3 at all redshifts. Our sur- 
vey mass limit typically lies above M* , thus our results 
are largely insensitive to the value of a. The HI mass func- 
tion model must specify the redshift-evolution of M*(z) and 
P(z). 

By integrating the mass function we can deduce the 
overall mass density of neutral hydrogen at a given redshift: 



pm(z) = <j>*{z)M*{z)T{a + 2), 



(11) 



where P denotes the Gamma function. We used the HIPASS 
survey mass function measurements to fix the value of 
Phi(0). In order to model the evolution of pm with red- 
shift, we followed AR05 and employed the measurements of 
the neutral gas density i n damped-Lyc t (DLA ) systems as 
a function of redshift by iPeroux et all (|200ll ). Specifically, 
we used the following functional form produced by AR05 
by fitting to DLA data points, making a somewhat ad-hoc 
correction for selection effects in the damped-Lya studies: 



fi H i(«) = N [1.813 - 1.473 (1 + z)- 2 - 31 ] 



(12) 



7 We n ote that due to missing fact ors of h in a plotting routine, 
Fig. 4 o f lAbdalla fc Rawlingsl 120051) is incorrect, and it should be 
replaced by FigfTJof this paper. 



where Qm(z) = pm{z) / p c {z = 0) is the HI mass density rel- 
ative to the critical density of the Universe (p c (z = 0) = 
2.7755 ft 2 10 11 M Q Mpc -3 ). The normalization constant is 
fixed by the zero-redshift HIPASS results. 

If we use Eq |12l to determine the redshift- variation of 
the normalization of the mass function, then the full evolu- 
tion model may be specified by the behaviour of the break in 
the Schechter function, M*(z). AR05 consider four different 
possibilities that, subject to the observational uncertainties 
in HI evolution detailed in the caption to Fig[2] span the 
plausible range of behaviours. We summarize these here, re- 
ferring the reader to AR05 for a more detailed discussion: 

• No-evolution model: (f>* and M* remain equal to the 
zero-redshift values measured by the HIPASS survey. 

• Model A: <j>* scales with z using the DLA results 
(Eq|T2]l with M*(z) = M*(0) = assumed constant. 

• Model B: cf>* scales with z using the DLA results. 
M*(z) = M*(0) D(z) 3 , where D(z) is the linear growth fac- 
tor at redshift z normalized such that D(0) = 1. 

• Model C: (j>* scales with z using the DLA results. The 
break in the mass function is additionally controlled by the 
cosmic star formation rate. 

In model C the break in the mass function evolves as 

M*(z) _ ( [atars(0)/O H l(0)] + 2 



M* (0) 



[fi s tar S (z)/OHl(z)]+2 



D{z) 



(13) 



where the fractional density in stars fi star s(.z) = 
p s ta,rs(z) / p c (z = 0) is deduc ed using the fits to the cos- 



mic st ar-formation history in IChoudhurv fc Padmanabhanl 
|2002l ) and there is an implicit assumption that ilm{z) = 
Qb\ 2 (z). Specifically, 



^sta 



Pc(z = 0) 



SFR(z') 



H{z')(l + z') 



Az (14) 



where H(z) is the Hubble parameter at redshift z and we 
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Figure 2. Measurements of the evolution of the HI mass den- 
sity relative to the critical de nsity of the Univer se, ttm(z), de- 
rived from HI emission studies l lZwaan et al,ll2005|) and studies of 
damped-Lyq absorption lines in distant quasars jProchaska et all 
l2005l ; lPeroux et al.ll2005h . We have corrected all data to our fidu- 
cial cosmological model. The solid line shows the evolutionary 
behaviour assumed in this paper, we have assumed following AR. 
This line adopts the maximal correction (factor 2) inferred from 
radio-selected samples for dust biasing against the optica l detec- 
tion of Hi-rich systems in absorption ijEllison et al.ll200lh . A re- 
cent study combining a bsorption-system results from optical- and 
radio-selected samples l|Pontzen fc Pettini 2008) suggests that the 
correction factor is unlikely to be much more than 1.3 subject to 
the caveat that a complete understanding of this factor requires 
more careful consideration of the effects of gravitational lensing, 
as well as dust, bias. 



take 



SFR(z) = 



0.13 



l + 6exp(-2.2,z) 



2 

2J5 



Xfld(z) 



x(z) 



(15) 



where x(z) is the co-moving distance to redshift z and 
Xfid(z) is the value of x(z) assuming a fiducial cosmology 
with Q m = 1 to convert from the cosmology assumed by 
IChoudhurv fc Padmanabhanl (|2002l ). The units of SFR(z) 
in Eq |15l are /i 2 M© yr _1 Mpc -3 . There is an implicit as- 
sumption that Qm(z) = Sln 2 { z )', although Qhi ~ 4S7jj 2 is 
more appropriate at low redshifts (Obreshkow & Rawlings 
2008), model predictions of the cosmic evolution in the de- 
gree of molecularization in galaxies (Obreshkow et al., in 
prep) suggest that a rough equality becomes a reasonable 
assumption at z ^ 1. Our default assumption will be evolu- 
tion Model C, following the recommendation of AR05. 



3.3 Survey parameters 

We consider HI emission-line surveys defined by four param- 
eters: 

• The telescope sensitivity as a fraction of that achieved 
by the SKA (the parameter / in Eq[2J). 

• The survey detection threshold n a in EqU 

• The telescope field-of-view at the HI emission-line fre- 
quency at zero redshift, FOV. 

• The assumed HI evolution model. 



The scaling of the observed wavelength of the HI emis- 
sion line with redshift can produce a corresponding scaling 
of the field-of-view FOV of the radio telescope, depending 
on the telescope design. We assume for now parabolic dishes 
with single-pixel feeds, such that FOV oc (l + z) 2 . This scal- 
ing may be used in a survey tiling scheme to produce an ef- 
fective increase with redshift in the integration time t of the 
HI survey, and this is incorporated into our model following 
AR05. 

We assume an overall survey duration equal to 1 yr 
covering approximately one hemisphere, 20,000 deg 2 . We do 
not consider trade-offs between total survey area and red- 
shift depth in this paper. 

The bandwidth of the radio telescope may not be suffi- 
cient to cover the wavelength range of the HI emission line 
up to the maximum observable redshift z max , necessitating 
repeat observations of each sky area. We use a parameter 
P to represent the ratio of the telescope bandwidth to the 
required frequency range of HI observations if we assume 
constant exposure time in each frequency block (in prac- 
tice more distant redshift slices will require longer integra- 
tions and an 'effective' value of /3 will result). For example 
if (3 x FOV = 10 deg 2 then a 20,000 deg 2 survey requires 
2000 pointings; if j3 = 1 each is of 4-hours duration, and 
this requires 1-year exposure time (adopting a 10 per cent 
overhead for calibration and other operational tasks). We 
note that, for a given survey duration and total solid an- 
gle, Eq[2] shows that the survey is identical if the combina- 
tion (P FOV f 2 /n a ) remains unchanged. For example, halv- 
ing the detection threshold n a is equivalent to doubling the 
field-of-view FOV or multiplying / by y/2. 

In subsequent calculations we adopt / = 0.15 as indica- 
tive of the sensitivity of an instrument representing 'phase-F 
of the SKA (which may be available in 2015). Although it is 
anticipated that such an instrument will have ~ 10 per cent 
of the final collecting area of the SKA, it is likely for prac- 
tical reasons to have a much larger fraction of its baselines 
within the compact core than is predicted by the baseline 
distribution function of Sec l3.ll Hence the sensitivity will be 
higher for fixed / as a greater fraction of baselines are useful 
for source detection. 



3.4 Clustering of Hi-emitting galaxies 

We assume that the Hi-emission galaxy power spectrum on 
large scales Pni(fc) is related to the underlying linear dark 
matter power spectrum Pii n (fc) by a constant bias factor b, 
as is observed for all classes of galaxy on large scales: 



Pni(k) = 6 2 P lin (fc) 



(16) 



Hi-emitting galaxies are known to be one of the least 
biased tracers of mass over-densities (i.e. to possess one 
of the lowest values of b) owing to their avoidance of 
dense regions such as galaxy clusters. The clustering of 
the zero-redshi ft HIP ASS survey has been measured by 
iBasilakos et al.l (|2007f ) by fitting a two-point correlation 
function parameterized by a power-law £(r) = (r/ro) 
with clustering length ro and slope 7. The results were ro = 
4.1 ± 0.3 h' 1 Mpc (in redshift space) and r = 3.1 ± 0.3 h' 1 
Mpc (in real space), using the fixed canonical value 7 = 1.8. 
This corresponds to a linear bias parameter b = 0.63. 

No equivalent measurement exists at higher redshift, 
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Figure 3. The dependence of the HI galaxy redshift distribu- 
tion dN/dz on the fractional sensitivity / of the radio telescope 
(where / = 1 for the SKA), assuming 1 yr of observing and 
(3 X FOV = 10 deg 2 . Other important assumptions are: (i) FOV 
scales with redshift as (1 + z) 2 ; (ii) the baseline distribution and 
galaxy size assumptions of Scc l3.1l arc adopted; (iii) HI evolution 
follows AR05 model C; (iv) only galaxies detected with a signifi- 
cance exceeding 10cr are considered. 



but we note that many classes of galaxy evolve in their clus- 
tering such that ro(z) « co nstant where rp is measured in 
co-moving co-ordinates (e.g. lLahav et al1l2002h . This forces 
a corresponding evolution in the linear bias b oc 1/D(z). 
Bias a ppears, however , to be a strong function of HI mass, 
in fact iBasilakos et~aL I (|2007l ) suggest that the more massive 
galaxies have a bias similar to optical galaxies, i.e. 6=1. 
For the purposes of this study and because of the lack of 
high-redshift measurements, we assume 6 = 1 independent 
of redshift; this is consistent with the clustering strength of 
objects near the break in the z — HI mass function. 

The clustering amplitude is significant in fixing the rel- 
ative contribution of cosmic variance and shot noise to the 
error in a measurement of the power spectrum. These two 
sources of error are equal (at a given Fourier scale k) when 
the quantity n x Pm(k) = 1, where n is the galaxy comoving 
number density. 



3.5 Redshift distributions 

Figs[3] and [4] plot the HI galaxy redshift distribution dN/dz 
for various different survey parameters assuming detection 
threshold n a = 10. The source density corresponding to 
nPm(k) = 1 (for our assumed HI clustering properties) 
is also indicated, using the characteristic scale k = 0.15 h 
Mpc - . The redshift at which dN/dz crosses this locus in- 
dicates the maximum redshift at which the error in the clus- 
tering measurements is dominated by cosmic variance rather 
than by shot noise. 

In Fig|3]we assume (3 FOV = 10 deg 2 , n CT = 10 and the 
preferred evolution model C, and plot dN/dz for fractional 
SKA areas / = 0.01, / = 0.15 and / = 1. In Fig g] we 
fix / = 1 and plot dN/dz for (3 FOV = 1, 10 and 100 deg 2 . 
Note that an identical redshift distribution is produced if the 
survey parameters are changed such that the combination 
{(3 FOV f fn a ) is left invariant. 



p FOV = 1 deg^ 
FOV = 10 deg 2 

FOV = 100 deg 2 



- 




Shot noise = cosmic variance 



2.5 



Redshift 



Figure 4. The dependence of the HI galaxy redshift distribution 
dN/dz on the quantity j3 X FOV, assuming 1 yr of observing and 
a telescope sensitivity / = 1. See the caption to Fig[3]for other 
assumptions. 



3.6 Effective survey volume 

A useful 'figure-of-merit' for the accuracy with which a 
galaxy survey spanning a range of redshifts can measure the 
power spectrum at a given scale k is the effective volume. 
This is the cosmic volume weighted by a redshift-dependent 
shot-noise factor: 



n(z)Pm(k, z) 



l + n(z)Pm{k,z) 



In our assessment we take k = 0.15 h Mpc -1 as the relevant 
scale, because this defines the approximate extent of the 
linear regime at z ~ 0.5 and also characterizes the regime of 
interest for the baryon oscillations. 

In Figs [5] and [6] we display the variation of the HI sur- 
vey effective volume with the telescope sensitivity / over 
the range from / = 0.01 to / = 1. For comparison we 
also indicate the approximate volume that is likely to be 
probed by ground-based surveys in optical wavebands on a 
timescale of 2015. We base this benchmark - V e s = 5 h~ 3 
Gpc 3 - on a projection of results from the Baryon Oscil- 
lation Spectroscopic Survey (BOSS) or the proposed Wide 
Field Multi-Object Spectrograph (WFMOS). This volume 
is equivalent to a survey of 10,000 deg 2 covering z < 0.7, 
or 2000 deg 2 spanning 0.5 < z < 1.5 (if we neglect the shot 
noise factor in Eq |17p . We note that assuming near-infra red 
spectroscopy, space-based surveys such as Euclid or JDEM 
would cover galaxies a redshift range of roughly z < 2.0 
(|Cimatti fc the SPACE collaborationll2008l ). hence a similar 
volume achievable by a full SKA. 

In Fig[5]we assume evolution model C and consider six 
different combinations of field-of-view and detection thresh- 
old. We note that the minimum sensitivity required by a 
radio telescope to equal the effective volume of our future 
optical survey benchmark (assuming n a = 10) are: 



(3 FOV = ldeg 2 -> / = 0.4 
(3 FOV = 10 deg 2 -> / : 
(3 FOV = 100 deg 2 / 



0.13 
= 0.04 
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This is consistent with the quantity FOV x f 2 characterizing 
the survey capability of radio telescopes. 

In Fig[S]we assume [3 FOV = 10 deg 2 and n a = 10 and 
display the dependence of the effective volume on the HI 
evolution model. The above conclusions regarding required 
collecting area are not strongly changed if a different evo- 
lution model is considered. In the most pessimistic case, an 
extra « 50% of sensitivity is required. 



3.7 Maximum survey redshift 

We used the criterion nPui(k) = 1 to estimate the max- 
imum redshift z max reachable in our model survey. Figs [7] 
and [S] display the dependence of this maximum redshift on 
the telescope collecting area in the same style as Figs[5]and 
[6] The HI radio surveys which match the performance of 
our benchmark future ground-based optical survey corre- 
spond to all-hemisphere surveys up to a maximum redshift 
z w 0.6 for an / = 0.15 radio array like the Phase-I SKA 
instrument. Surveys with a full SKA would access redshifts 
significantly beyond z = 1, perhaps reaching z ~ 2. Again, 
these conclusions do not depend strongly on the assumed HI 
evolution model. 



4 SIMULATED DARK ENERGY 

MEASUREMENTS USING BARYON 
OSCILLATION FITTING 

4.1 Survey power spectra 

We divided each model HI survey into redshift slices of width 
Az — 0.2 up to the maximum redshift z max defined in Sec l3.7l 
(rounding up to the next multiple of Az). For each redshift 
slice we modelled the error in the measurement of the galaxy 
power spectrum Phi (k) in each of a series of Fourier bins of 
width Afc = 0.01 h Mpc" 1 : 



SP m (k,z) = — 



1 + 



n(z) P m (k,z) 



(18) 



where n(z) is the HI galaxy number density at redshift z, 
and m is the number of Fourier modes appearing in the bin 
of width Ak: 



1 



V 



m = — X -. rrr 

2 (2tt) 3 



x Ank Ak 



(19) 



where V is the total survey volume. The power spectrum 
errors in each of the redshift slices are displayed in the panels 
of Fig|9] for four example surveys: 



• Parameter set 1: / = 

• Parameter set 2: / 



10. 



• Parameter set 3: / 



10. 



• Parameter set 4: / 



1, (3 FOV = 10 
0.15, I3FOV = 

0.01, (3FOV = 

0.15, (3FOV = 



10 deg 
10 deg 2 
10 deg 2 



n a = 10. 

2 „ 



In each case the model power spectra and errors have 
been divided by a 'reference' power spectrum which pos- 
sesses the same broadband shape and amplitude as the 
model power spectrum but has the BAOs erased, i.e. the 



'no- wiggles' spectrum of lEisenstein fc Hul (|l998l ). Present- 
ing the data in this manner illustrates the approximate sig- 
nificance with which the baryon oscillations are delineated 
by these surveys. The different redshift slices are offset ver- 
tically for clarity, and data points are plotted up to a max- 
imum value k = 0.2 h Mpc -1 , which is an estimate of the 
maximum extent of the linear regime at low redshifts, and 
which we then conservatively approximate as being inde- 
pendent of redshift. We note that an / = 0.15 SKA could 
produce a high-fidelity measurement of baryon oscillations 
in each of a series of redshift slices up to z « 0.8, if equipped 
with a field-of-view (3 FOV > 10 deg 2 . A full SKA could re- 
turn exquisite measurements of the galaxy power spectrum 
up to z w 2. 

4.2 Fitting formula results 

The BAOs in the galaxy power spectra act as a cosmological 
standard ruler which enables the cosmic distance and expan- 
sion rate at ea ch survey redshift slice to be measured with 
high precision ([Blake fc Glazebroo k 2003; S ep fc Eisensteinl 
|2003| ; Ih"u fc Haimanll2003T ). Roughly speaking, oscillations in 
the tangential direction of the dataset can be used to accu- 
rately map the cosmic distance to each redshift slice. Radial 
osci llations return the H ubble parameter H(z) at each slice. 

iBlake et all (|2006h devised fitting formulae which al- 
low the survey parameters - central redshift z, volume V 
and galaxy number density n - to be converted into ap- 
proximate accuracies of measurement of the tangential and 
radial baryon oscillation scales . These formulae are based 
on t he procedure discuss e d by IBlake fc Glazebrook! (|2003l ) 
and iGlazebrook fc Blakei l|2005h . which is designed to pro- 
duce robust an d conservative est i mates of the standard ruler 
accuracies fsee lSeo fc Eisensteinl (|2007l ) for a more recent fit- 
ting formula). 

We used the fitting formulae of IBlake" et all (|2006h to 
determine the tangential and radial standard ruler accura- 
cies in each of a series of redshift slices of width Az = 0.2 
for a variety of surveys varying collecting area /, detection 
threshold n CT and field-of-view (3 FOV. We assume HI evo- 
lution model C. The results are plotted in Fig llOl The solid 
and dashed lines plot the tangential and radial accuracy, re- 
spectively, which are approximately equal to the accuracy 
with which the quantities x(z) and H(z) may be inferred 
in each redshift slice. This plot serves to emphasize the im- 
portance of a large field-of-view in extending the redshift 
range over which an HI survey can map cosmic distances. 
Future radio telescopes can produce per-cent-level measure- 
ments of x(z) and H(z) over a range of redshifts, provided 
that (3 FOV f 2 > 1. 

These measurements of cosmic distances and expansion 
rates may be readily converted into determinations of the 
parameters of the dark energy model. This is illustrated by 
Figs llll and ll2l corresponding to two different parameteriza- 
tions of the dark energy equation of state: (1) w(z) = ui C ons 
(Fig[TT]| and (2) w(z) = wo + w a (l — a), where a = l/(l + z) 
fFig |12|l . We assume Gaussian priors on the other relevant 
cosmological parameters - fi m and h - parameterized as 
a(Qmh ) = 0.004 and cr(fi m ) = 0.01. This is representative 
of the combination with measurements of the Cosmic Mi- 
crowave Background anisotropies by the forthcoming Planck 
satellite. In detail, we assume that the baryon oscillations 
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Fraction of SKA sensitivity 

Figure 5. The dependence of the effective HI survey volume achieved in 1 yr of observing on the fraction of SKA sensitivity, for various 
different values of field-of-view (/3 X FOV) and detection threshold n a . See the caption to Fig|3]for other assumptions. A rough benchmark 
is plotted for future ground-based redshift surveys in optical and near-infrared wavebands. 




Fraction of SKA sensitivity 

Figure 6. The dependence of the effective HI survey volume achieved in 1 yr of observing on the fraction of SKA sensitivity, for different 
HI evolution models. We assume j3 X FOV = 10 deg 2 ; see the caption to Fig[3]for other assumptions. A rough benchmark is plotted for 
future ground-based redshift surveys in optical and near-infrared wavebands. 
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Fraction of SKA sensitivity 

Figure 7. The dependence of the maximum HI survey redshift z max achieved in 1 yr of observing on the value of /3 X FOV and detection 
threshold n CT . See the caption to Fig|3]for other assumptions. 
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Fraction of SKA sensitivity 

Figure 8. The dependence of the maximum HI survey redshift z max achieved in 1 yr of observing on the fraction of SKA sensitivity, for 
different HI evolution models. A telescope with /3 X FOV = 10 deg 2 is assumed; see the caption to Fig[3]for other assumptions. 
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Figure 9. Power spectrum simulations for a series of redshift slices for four different HI survey parameter sets. The power spectra are 
divided by a smooth reference spectrum to emphasize the baryon acoustic oscillations. Power spectra for different redshift slices are offset 
vertically for clarity. Points ar e plotted up to a maximum value k = 0.2 h Mpc - 1 which is an approximate estimate of the extent of the 
linear regime at low redshifts l|Blake fc Glazebrooldl2003h . 
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Figure 10. Simulated standard ruler accuracies in redshift slices of width Az = 0.2 for a series of survey parameter sets defined by 
fractional SKA area /, threshold for source detection n CT , and field-of-view j3 X FOV. The solid lines indicate the simulated accuracy of 
measurement of the tangential acoustic scale; the dashed lines display the accuracy of the radial scale. 
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Figure 13. We plot here the accuracies in determining the an- 
gular diameter distance and the Hubble parameter from some 
model optical surveys, using the fitting formula described and wc 
over-plot results for radio surveys. We have kept the binning in 
dz=0.2 redshift slices The optical surveys we use are: 0.2 j z | 0.4, 
10000 deg 2 , 0.4 i z i 0.6, 10000 deg 2 , 0.6 j z i 0.8, 2000 deg 2 , 0.8 
j z j 1.0, 2000 deg 2 and 1.0 j z j 1.2, 2000 deg 2 . nP = 2 for all 
surveys. nP = 2 refers to the product of the galaxy number den- 
sity n in the redshift slice and the amplitude of the galaxy power 
spectrum P at k=0.2 (i.e. it is independent of galaxy bias) and 
ensures shot noise is significantly less than cosmic variance. The 
comparison radio surveys all have /3 FOV = 10 deg 2 , duration of 
one year, area of 20,000 deg 2 . The solid lines (solid circles) are 
angular diameter distances accuracies and the dashed lines (open 
circles) are Hubble parameter accuracies. 



provide a measure of x(z)/s and H(z)~ 1 / s, where s is the 
sound horizon at recombination, and fold the uncertainties 
in Q m and h into the error in s. 

An / = 0.15 SKA can produce a strong measurement 
of a constant dark energy equation of state (o-(ui CO ns) < 0.1) 
if (3 FOV > 10 deg 2 . A full SKA produces <r(w coa3 ) w 0.025 
for (3 FOV — 10 deg 2 . For an evolving equation of state 
parameterized by (wo,w a ), measurements at higher redshift 
assume increased importance. Here, an / = 0.15 SKA can 
produce a good measurement of the evolving term (a(w a ) < 
1) if (3 FOV > 50 deg 2 assuming n a = 10. A more aggressive 
threshold n CT = 5 doubles the accuracy of the dark energy 
measurement in this case. A full SKA with f3 FOV = 10 deg 2 
allows parameter measurements of superb quality: <j(wo) fs 
0.05 and a(w a ) « 0.2. In Sec[S]we will compare these results 
with a full MCMC calculation. 

We also plot in Fig |13l the accuracies in determining the 
angular diameter distance and the Hubble parameter from 
some model optical surveys using the fitting formula as we 
already described and over plotted onto results for radio 
surveys. This plot gives us an idea of the relative merit of 
optical surveys compared to radio surveys for cosmological 
parameter estimation. A binning of 0.2 in redshift was used 
and a number density for optical surveys was assumed to 
be nP — 2 (see Eqn ll8[) . the product of the galaxy num- 
ber density n in the redshift slice and the amplitude of the 
galaxy power spectrum P at k — 0.2 (i.e. it is independent of 
galaxy bias) - nP — 2 ensures shot noise is less than cosmic 
variance. 



5 SIMULATED DARK ENERGY 

MEASUREMENTS USING A FULL POWER 
SPECTRUM ANALYSIS 

We now perform simulations using the full power spectrum 
information to measure the dark energy model. The broad- 
band shape of the power spectrum depends mainly on the 
quantity Q m h, and redshift-space distortions contain infor- 
mation on the large-scale bias of the galaxy population. 
Hence the power spectrum shape helps to break degenera- 
cies in the cosmological model and sha rpen our dark energy 
measu rements. We refer the reader to lAbdalla fc Rawlingsl 
(|2007l ) for a full description of our MCMC forecast method. 

In this Section we assume a survey performed by an 
SKA with (3 x FOV = 10 deg 2 . Therefore, a one-year sur- 
vey corresponds to 4 hours integration time per pointing. 
There is a straight-forward trade-off between collecting area 
and integration time; when comparing different collecting 
areas we specify for reference the integration time required 
to achieve an equivalent depth with the full SKA (/ = 1). 

5.1 Effects of survey depth on parameter forecasts 

We present in this subsection the effect that different col- 
lecting areas (or equivalently different integration times with 
a full SKA) have on the accuracy of cosmological param- 
eter measurement for a 20,000 deg 2 survey. Initially we 
base the forecasts only on the SKA large-scale structure 
data and take a standard ACDM framework with an extra 
parameter w, a constant equation-of-state of dark energy. 
Our fiducial model is: {Sl6,f7 c ,i«,/i,n s ,(Tg,r} = {0.04,0.26,- 
1.0,0.72,1.0,0.9,0.09}. The fiducial bias of each redshift bin 
is b = 1, and we use three redshift bins in the analysis. We 
have tested that for a constant-iu model there is negligi- 
ble benefit from analyzing more redshift bins, although we 
note that this may not be the case for a more complex dark 
energy parameterization. We present in Table [T] our fore- 
cast errors in cosmological parameters recovered from the 
MCMC chains. 

Using the large-scale structure data alone, our measure- 
ment error in w degrades sharply at a collecting area below 
/ ~ 0.35 (or a 30-min exposure using the full SKA). This is 
because below this approximate threshold we are no longer 
able to break the degeneracies between the fitted parameters 
and we require additional data such as the CMB. 

5.2 A combined experiment 

We next forecast the cosmological parameter measurements 
achieved by combining the SKA large-scale structure sur- 
vey with the CMB ( Planc k) dataset using the method of 
lAbdalla fc Rawlingsl (|2007h . We examine how an / = 0.15 
SKA survey and an / = 1 SKA survey would benefit from 
the inclusion of CMB data. In order to combine the two 
datasets we add the log(likelihood) distributions assuming 
the datasets are uncorrelated. 

We plot in Fig |15l the parameter measurements found 
for the combination of full SKA survey and CMB data. 
There are no significant degeneracies in any of the parameter 
pairs, as opposed to (for example) the w — erg plane plotted 
in Fig |14l from the SKA survey alone (see also Fig |16[) . The 
CMB provides an accurate amplitude of fluctuations at the 
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Figure 11. Simulated dark energy measurements of a constant cquation-of-state w(z) = w conB for a series of surveys defined by fractional 
SKA sensitivity /, threshold for source detection n CT , and field-of-view f3 X FOV. 1-cr confidence ellipses are plotted in the parameter 
space of (n m , tu cons ). Marginalized errors are shown for the parameter ui cons . The analysis uses the measurements of x(z) and H (z) using 
tangential and radial baryon oscillations, as displayed in Fig llOl together with Gaussian priors a(Q m h 2 ) = 0.004 and cr(n m ) = 0.01, 
representative of combining with Planck satellite observations of the Cosmic Microwave Background. 




Figure 12. Simulated dark energy measurements of a redshift-dependent equation-of-state w(z) = too + tl)«(l — a) for a series of surveys 
defined by fractional SKA sensitivity /, threshold for source detection n CT , and field-of-view (3 X FOV. 1-cr confidence ellipses are plotted 
in the parameter space of (wo, w a ). Marginalized errors are shown for wo and w a . The analysis uses the measurements of x(z) and H (z) 
using tangential and radial baryon oscillations, as displayed in Fig llOl together with Gaussian priors cr(n m /i 2 ) = 0.004 and cr(Q m ) = 0.01, 
representative of combining with Planck satellite observations of the Cosmic Microwave Background. 



© 0000 RAS, MNRAS 000, 000-000 



14 F. B. Abdalla et al. 





\ 




















/ 








X 


• 




(§) 











4 0.0410.042 0.25 0.255„ 0.26 0.2650.715 Q.72 0.725 0.99 



1.02 0.895 0.9 0.905 




0.995 1 1.005 1.01 0.8B 0.9 0.92 0.94 0.96 0.16 0.1^ 



0.262 

0.261 

; u 0.26 
0.253 

2:>, 



0.722 
0.721 
c 0.72 
0.71 
0.718 
0.717 



0.004 
002 



0.175 
0.17 
" 0.165 
0.16 

0.155 

































® 




(®) 


• 








^ Xv ~x^x 








* 








.0) 





0.039R03980. 040.04020404 0.258 



0^6 0.262 0.718 0.72 0.7 



1.005 0.896 0.898 0.9 0.902 0.903.1550.160.1650.170.175 



Figure 15. Analysis of the accuracy of a combination of future SKA (red), CMB (blue) and combined (green) experiments in measuring 
cosmological parameters. The contours represent one and two sigma confidence regions. The combination of SKA and CMB experiments 
will measure all current cosmological parameters to below the one per cent level, including the dark energy parameter w. The full SKA 
will produce a measurement of w with an accuracy close to 0.5%. 
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Table 1. The forecast errors in cosmological parameters as a function of survey depth. The top row indicates the telescope sensitivity 
(for a 4-hr exposure); the second row lists the equivalent exposure time with a full SKA (/ = 1). The third row details the approximate 
maximum redshift achieved, i.e. where nP = 1. For small integration times the errors in the dark energy equation-of-state w are large 
because the baryon oscillations are not measured accurately. 
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Figure 14. The posterior in the (erg, w) plane for full SKA sur- 
veys with different exposure times. The forecasts improve signif- 
icantly from 6 mins exposure time (equivalent to 4 hrs with an 
/ = 0.15 SKA) to 4 hrs exposure time with the full SKA. When 
we reach depths equivalent to a 1-hr exposure time with the full 
SKA, there are sufficient galaxies per unit volume and redshift to 
break the degeneracies and produce an accurate measurement of 
w. 




Figure 16. We illustrate one of the main degeneracies resolved 
when an HI survey dataset obtained with the / = 0.15 SKA is 
combined with a CMB (Planck) dataset. The degeneracy between 
as ancl io is broken and allows for a very good measurement of 
the dark energy parameter, even with a relatively shallow survey 
covering 20,000 deg 2 . The samples are for independent MCMC 
chains for an / = 0.15 SKA and for Planck. They are colour-coded 
according to the value of h, as shown in the side-bar. 



last scattering surface, together with a measurement of the 
angular diameter distance to this surface, adding comple- 
mentary information to the SKA dataset. In Table[5]we can 
see that there is a significant improvement in the parameter 
errors compared to those found in Table [T] The accuracy of 
measurement of the w parameter is below 1 per cent. 

We also conclude from the simulations that most of the 
constraining power with regard to the parameter n s origi- 
nates from the CMB, and the LSS dataset adds little. This 
is because a high-resolution CMB experiment with a large 
multipole range (such as Planck) will measure CMB fluc- 
tuations across a large range in spatial scales. This picture 
could change if one were able to model non-linearities in the 
large-scale structure including a reliable biasing model at 
small scales. We surmise that the CMB data will be very 
useful in constraining some parameters, e.g. n a , whilst LSS 
surveys will play an increasing role in the determination of 
other parameters, such as as and w (although the capacity 



to measure as rests on a good model for the galaxy bias 
from redshift-space distortions). 

We plot in Fig |16l the MCMC samples obtained for an 
/ = 0.15 'phase I' SKA combined with the future CMB data. 
In the absence of systematic errors this would yield errors 
in w of approximately 1 per cent. This forecast worsens by a 
factor of 2 if we assume that redshift-space distortions can- 
not be used to determine the galaxy bias factor, as discussed 
in the next Section. 



5.3 Normalization factors 

The results in Sec[5] rest on the assumption that redshift- 
space distortions can constrain the galaxy bias factor via 
large -scale bulk flows (|Cuzzo et al. 1 120081 ; iPercival fc White] 
l200cl ). This information can also be provided by other 
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Figure 17. Wc demonstrate how the accuracy in measuring the 
value of the dark energy parameter w worsens if we cannot re- 
cover the galaxy bias information. In this case the distribution 
of MCMC chains changes from the black points to the coloured 
points, a degeneracy between erg and w opens up, and the error 
in w degrades from 1% to 2.5%. The MCMC samples without 
bias information are colour-coded according to the value of h, as 
shown in the side-bar. 



probes which we have not used here. Firstly, the analy- 
sis of the small-scale clustering in the framework of the 
galaxy halo model makes predictions for the large-scale 
bias and redshift-space effects. Secondly, higher-order clus- 
tering statistics y ield information on the bias of galaxies 
I Verde et alj|2002l ). Reasonably accurate knowledge of the 
bias parameter(s) may be gleaned from these tools. The 
scale-dependence of bias may be tested by selecting galax- 
ies with different c lustering propert ies located in the same 
comoving volume (|Cole et alj l2006h , with the aim of de- 
coupling the cosmological signal from the galaxy bias. 

On the other hand, we should be cautious in making 
these assumptions, since it is not yet clear that the desired 
level of systematic error control can be achieved. In order to 
investigate the importance of the additional bias information 
we performed a simulation considering only the isotropized 
power spectrum obtained with an / = 0.15 SKA survey in 
conjunction with the Planck CMB data, thereby removing 
the information associated with the redshift-space distor- 
tions. 

We find that the parameter which suffers most from the 
loss of the bias information is erg , which is unsurprising given 
that it is directly related to the height of the power spec- 
trum. In Fie ll7l we plot results from simulations, both with 
and without bias measurements from redshift-space distor- 
tions. We find that without the bias information the accu- 
racy of measurement of w is reduced from the one per cent 
level to roughly the two and a half per cent level, and there 
remains a residual a$-w correlation. 



Parameter / = 0.15 SKA (1) / = 0.15 SKA (2) / = 1 SKA 



Q b 0.00024 0.00044 0.00016 

U c 0.0016 0.0021 0.00086 

w 0.010 0.022 0.0062 

n a 0.0024 0.0028 0.0021 

o-g 0.0026 0.010 0.0013 

h 0.0019 0.0035 0.0010 

t 0.0042 0.0041 0.0036 



Table 2. The forecast errors in cosmological parameters from 
simulations of an SKA HI survey plus CMB (Planck) data. The 
first column of results is for a 'phase-F (/ = 0.15) SKA and cor- 
responds to a full analysis including the redshift-space distortions 
effect. The second column assumes that this information is un- 
reliable and marginalizes over the large-scale anisotropic power. 
The final column compares these results to a full SKA survey 
including redshift-space distortions. 

5.4 Comparison with BAO-only results 

We can compare these MCMC calculations with our more 
approximate results presented in Sec|4] In the approxima- 
tion of SecU we only utilize the distance information en- 
coded in the BAOs, not the full power spectrum shape. In 
addition we incorporate the CMB data only as 'representa- 
tive' priors in the values of Q m h 2 and f2 m rather than by 
calculating the full parameter degeneracies, and we do not 
use information encoded in the growth of structure. 

In Sec(3]we forecast an error ti(w) — 0.026 for an HI 
survey with the full SKA with j3 FOV = 10 deg 2 . The error 
determination in Table [2] is ~ 4 times smaller. We suggest 
that approximately half of this improved performance is the 
result of utilizing the growth factor measurements, and the 
remainder is the result of using the shape of the power spec- 
trum to measure the values of f2 m and h more accurately 
hence sharpening our capacity to utilize the BAO distance 
information. 



6 REMARKS ON SYSTEMATIC ERRORS IN 
FUTURE RADIO SURVEYS 

The forecasts in this paper assume that the HI galaxy clus- 
tering measurements are not compromised by systematic 
errors from spurious sources or artifacts in the data. We 
present results for 10-a and 5-a source catalogues; however, 
there is a wealth of extra information contained in the lower- 
sigma peaks in our derived radio images, and our analysis 
is conservative. In principle it might be possible to analyze 
the raw intensity data cube without the intermediate step of 
extracting a catalogue and thereby greatly enhance the max- 
imum redshift and number density of our HI galaxy survey 
jChang et al.ll2008l ). 

A concern with this approach is our ability to quan- 
tify accurately the power spectrum of the noise. Given the 
dynamic range limitations of radio interferometry, the noise 
will be highly variable across the sky and depend on a host of 
subtle calibration issues; it is certainly likely to be enhanced 
in the vicinity of bright radio continuum sources imprinting 
characteristic scales in the measurement. The imprint of this 
variation will require very careful modelling. It is also possi- 
ble that, in contrast to CMB observations, the noise in the 
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HI data cube may be non-Gaussian. A detailed investigation 
of these issues is beyond the scope of this paper. 

Another source of systematic error in the galaxy power 
spectrum is the uncertainty in the galaxy bias factor. This 
can be investigated by measuring the clustering properties 
of galaxy sub-samples (split for example by HI mass) across 
which the imprint of cosmological parameters should be 
identical. This programme can be carried out effectively for 
HI surveys owing to the large number density of sources. 
This is illustrated by Fig |18l - the solid lines show the value 
of the quantity nP for different models, assuming a 1-yr 
SKA survey with j3 FOV ~ 10, where n is the galaxy num- 
ber density and P is the power spectrum amplitude. The 
product nP > 10 over a large redshift range, implying that 
we can measure a power spectrum without serious shot noise 
error for more than ~ 10 subsets of galaxies. An aggressive 
(few a) detection threshold may require more complicated 
treatments of galaxy bias. Low-threshold catalogues suffer 
heavily from Eddington bias effects that (considering for 
simplicity one fixed redshift) will yield a spuriously wide 
range of intrinsic HI mass (and hence galaxy bias) at the 
fixed detection HI line detection threshold. 

Because this survey is 'over-sampled' at low redshifts it 
is possible to reduce significantly the performance of the tele- 
scope without compromising much of the cosmological sci- 
ence. In Fig |18l we consider the case where the sensitivity de- 
creases with the square of the frequency above 800 MHz; the 
result is the dashed line and this could yield important de- 
sign considerations: for example a quadratic drop-off in the 
sensitivity of aperture arrays arises naturally if critical sam- 
pling of the wavefront (having antenna elements spaced at 
least every half wavelength) becomes relaxed below a wave- 
length, say that corresponding to 800 GHz, meaning that 
at these higher fre quencies the array becomes sp arse rather 
than close packed (jBraun fc van Capelle nl l2006l ). Gaining a 
clearer idea of the frequency ranges where noise can be com- 
promised depends on a better measurement of the HI mass 
function at z ~ 1, emphasizing the importance of the stack- 
ing experiments mentioned in Sec(2] 

In our baseline model the full SKA can survey the z < 2 
Universe on the timescale of a year, based on technologies 
such as a phased array or dishes with focal-plane phased ar- 
rays in order to achieve the required fields of view. At lower 
frequencies, sparsely-sampled dipoles may be the technology 
of choice. In this case the scaling of field-of-view with fre- 
quency can be set to the user requirements provided that the 
computing power is available to process the required num- 
ber of beams. Our assumed (1 + z) 2 scaling of the FOV is 
just one choice, that needs to be optimized with detailed 
simulations. At the lowest frequencies it may be desirable to 
have the physical collecting area rising as a sharp function 
of wavelength to ensure a mapping speed that is constant 
in redshift. For a sparse array this is naturally achieved via 
the A 2 variation in the collecting area of antennas. However 
it may still be necessary to build physically more antennas 
in this frequency range, or to use the processing system to 
shape how the FOV varies with frequency, and hence red- 
shift, by, for example, forming more beams at the lower fre- 
quencies to yield greater instantaneous sky coverage. This is 
subject to the fundamental limitation of the finite sky area 
observable at moderate zenith angle from an aperture array 
station. 
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Figure 18. The solid lines display the variation of the value of 
nP with redshift for the 3 different HI evolution models. For all 
models the value of nP remains above 1 for z < 1, implying 
that we can measure the galaxy power spectrum for many sub- 
samples and hence constrain the systematic effects of galaxy bias. 
The dashed lines illustrate the degradation which occurs when the 
telescope sensitivity worsens as the frequency squared in the high 
frequency range, i.e. at low redshift. 

7 CONCLUSIONS 

We have forecast the cosmological parameter measurements 
that could result from a combination of future Hi-emission 
galaxy surveys with the SKA and CMB (Planck) data. The 
main degeneracies currently present in cosmological data 
would be broken by such new surveys. Current parameters 
would be measured to below the 1% level and other param- 
eters such as neutrino masses will be accessible. 

We find that an HI survey with a telescope possess- 
ing / = 0.15 of the SKA sensitivity, together with a future 
CMB experiment, will be able to determine the dark energy 
equation of state with an accuracy close to one per cent. 
We have examined the origin of this information and found 
that the robust signature of 'geometry' (the projection of 
the BAOs) is producing significant power. The remaining 
information originates from measurements of the growth of 
cosmic structure that can be inferred by looking at redshift- 
space distortions at large scales and thereby breaking the 
degeneracy with galaxy bias. 

We caution that the evolutionary model considered here 
adopts a slightly stronger evolution in the cosmic mass den- 
sity of HI (see Fig. 2) than i s inferred from the most re- 
cent studies of HI absorption (|Pontzen fc Pettirul 120081 ) . so 
the observing times required may be underestimated by a 
small factor (assuming most of the number counts are com- 
ing from the exponential tail of the mass function near the 
break, which should be the case in the high redshift region; 
in order to have ~ 1.5 times more sources, one would need 
to detect sources ~ 1.7 times the mass, which would require 
~ 2.8 the exposure time roughly). We emphasise the impor- 
tance of stacking experiments in Sec(2] 

A survey with the full SKA would be able to produce 
a measurement of the dark energy parameter with an ac- 
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curacy close to half a per cent. In this case the addition of 
the CMB data does not produce such a substantial increase 
in performance. However, the SKA and CMB data remain 
complementary in the sense that the galaxy survey data is 
better able to measure the height of the power spectrum 
(hence ag) and the CMB data is most powerful for deter- 
mining the spectral index of scalar fluctuations (n s ). 

We conclude that a radio telescope with a sensitivity 
between / = 0.15 and / = 1 of the SKA would be a cosmo- 
logical tool capable of transformational science, particularly 
in the field of dark energy. 



ACKNOWLEDGMENTS 

We thank Kathryn Nicklin and Danail Obreschkow for prac- 
tical help with the work of this paper and for useful dis- 
cussions. FBA acknowledges support from the Leverhulme 
Foundation via an Early Careers Fellowship. This work has 
been undertaken as part of the Square Kilometre Array De- 
sign Study (SKADS) financed by the European Commission. 



REFERENCES 

Abdalla F. B., Rawlings S., 2005, MNRAS, 360, 27 
Abdalla F. B., Rawlings S., 2007, MNRAS, 381, 1313 
Basilakos S., Plionis M., Kovac K., Voglis N., 2007, MN- 
RAS, 378, 301 
Blake C, Glazebrook K., 2003, ApJ, 594, 665 
Blake C, Parkinson D., Bassett B., Glazebrook K., Kunz 

M., Nichol R. C, 2006, MNRAS, 365, 255 
Blake C. A., Abdalla F. B., Bridle S. L., Rawlings S., 2004, 

New Astronomy Review, 48, 1063 
Boughn S. P., Crittenden R. G, 2005, New Astronomy 

Review, 49, 75 
Braun R., van Capellen W., 2006 

Cabre A., Gaztanaga E., Manera M., Fosalba P., Castander 
F., 2006, mnras, 372, L23 

Catinella B., Haynes M. P., Giovanelli R., Gardner J. P., 
Connolly A. J., 2008, ApJ Lett., 685, L13 

Chang T.-C, Pen U.-L., Peterson J. B., McDonald P., 2008, 
Physical Review Letters, 100, 091303 

Choudhury T. R., Padmanabhan T., 2002, ApJ, 574, 59 

Cimatti A., the SPACE collaboration 2008, Experimental 
Astronomy, pp 37 — h 

Cole S., Sanchez A. G, Wilkins S., 2006, ArXiv Astro- 
physics e-prints 

Cruz M. J., Jarvis M. J., Rawlings S., Blundell K. M., 2007, 
mnras, 375, 1349 

Dalton G. B., Lewis I. J., Bonfield D. G, Holmes A. R., 
Brooks C. B., Lee H., Tosh I. A. J., Froud T. R., Patel 
M., Dipper N. A., Blackburn C, 2006, in Ground-based 
and Airborne Instrumentation for Astronomy. Edited by 
McLean, Ian S.; Iye, Masanori. Proceedings of the SPIE, 
Volume 6269, pp. 62694A (2006). Vol. 6269 of Presented at 
the Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference, The UK FMOS spectrograph 

Duffy A. R., Battye R. A., Davies R. D., Moss A., Wilkin- 
son P. N., 2008, MNRAS, 383, 150 

Dunkley J., Komatsu E., Nolta M. R., Spergel D. N., Lar- 
son D., Hinshaw G, Page L., Bennett C. L., Gold B., 



Jarosik N., Weiland J. L., Halpern M., Hill R. S., Kogut 
A., Limon M., Meyer S. S., Tucker G. S., Wollack E., 
Wright E. L., 2008, ArXiv e-prints, 803 

Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 

Eisenstein et al. 2005, ApJ, 633, 560 

Ellison S. L., Yan L., Hook I. M., Pettini M., Wall J. V., 
Shaver P., 2001, 379, 393 

Freedman W. L., Madore B. F., Gibson B. K., Ferrarese 
L., Kelson D. D., Sakai S., Mould J. R., Kennicutt R. C, 
Ford H. C, Graham J. A., Huchra J. P., Hughes S. M. G, 
Illingworth G. D., Macri L. M., Stetson P. B., 2001, ApJ, 
553, 47 

Giannantonio T., Crittenden R. G, Nichol R. C, Scranton 
R., Richards G. T., Myers A. D., Brunner R. J., Gray 

A. G, Connolly A. J., Schneider D. P., 2006, prd, 74, 
063520 

Glazebrook K., Blake C, 2005, ApJ, 631, 1 

Glazebrook K., Blake O, Couch W., Forbes D., Drinkwater 
M., Jurek R., Pimbblet K., Madore B., Martin C, Small 
T., Forster K., Colless M., Sharp R., Croom S., Woods 
D., Pracy M., Gilbank D, Yee H., Gladders M., 2007, in 
Metcalfe N., Shanks T., eds, Cosmic Frontiers Vol. 379 
of Astronomical Society of the Pacific Conference Series, 
The WiggleZ Project: AAOmega and Dark Energy, pp 
72-+ 

Guzzo et al. 2008, Nat, 451, 541 

Hill G. J., Gebhardt K., Komatsu E., Drory N., MacQueen 
P. J., Adams J., Blanc G. A., Koehler R., Rafal M., Roth 
M. M., Kelz A., Gronwall C, Ciardullo R., Schneider 

D. P., 2008, ArXiv e-prints, 806 

Hu W., Haiman Z., 2003, Phys. Rev. D, 68, 063004 
Ishibashi A., Wald R. M., 2006, Classical and Quantum 

Gravity, 23, 235 
Komatsu E., Dunkley J., Nolta M. R., Bennett C. L., Gold 

B. , Hinshaw G, Jarosik N., Larson D., Limon M., Page L., 
Spergel D. N., Halpern M., Hill R. S., Kogut A., Meyer 
S. S., Tucker G. S., Weiland J. L., Wollack E., Wright 

E. L., 2008, ArXiv e-prints, 803 
Lahav et al. 2002, MNRAS, 333, 961 

Mena O., Santiago J., Weller J., 2006, Physical Review 

Letters, 96, 041103 
Peroux C, Irwin M. J., McMahon R. G, Storrie-Lombardi 

L. J., 2001, Astrophysics and Space Science Supplement, 

277, 551 

Percival W. J., Nichol R. C, Eisenstein D. J., Frieman 
J. A., Fukugita M., Loveday J., Pope A. C, Schneider 
D. P., Szalay A. S., Tegmark M., Vogeley M. S., Weinberg 
D. H., Zehavi I., Bahcall N. A., Brinkmann J., Connolly 
A. J., Meiksin A., 2007, ApJ, 657, 645 

Percival W. J., Nichol R. C, Eisenstein D. J., Weinberg 
D. H., Fukugita M., Pope A. C, Schneider D. P., Szalay 
A. S., Vogeley M. S., Zehavi I., Bahcall N. A., Brinkmann 
J., Connolly A. J., Loveday J., Meiksin A., 2007, ApJ, 
657, 51 

Percival W. J., White M., 2008, ArXiv e-prints, 808 
Peroux C, Dessauges-Zavadsky M., D'Odorico S., Sun Kim 

T., McMahon R. G, 2005, MNRAS, 363, 479 
Pontzen A., Pettini M., 2008, ArXiv e-prints 
Prochaska J. X., Herbert-Fort S., Wolfe A. M., 2005, ApJ, 

635, 123 

Rasanen S., 2004, Journal of Cosmology and Astro-Particle 
Physics, 2, 3 



© 0000 RAS, MNRAS 000, 000-000 



Dark Energy with future HI surveys 

Rassat A., Land K., Lahav O., Abdalla F. B., 2006, ArXiv 

Astrophysics e-prints 
Rawlings S., Abdalla F. B., Bridle S. L., Blake C. A., Baugh 

C. M., Greenhill L. J., van der Hulst J. M., 2004, New 

Astronomy Review, 48, 1013 
Refregier A., the DUNE collaboration 2008, ArXiv e-prints 
Riess et al. 2006, ArXiv Astrophysics e-prints 
Seo H.-J., Eisenstein D. J., 2003, ApJ, 598, 720 
Seo H.-J., Eisenstein D. J., 2007, ApJ, 665, 14 
Simpson F., Bridle S., 2006, prd, 73, 083001 
Skordis C, Mota D. F., Ferreira P. G., Bcehm C, 2006, 

Physical Review Letters, 96, 011301 
Spergel et al. 2007, apjs, 170, 377 

Taylor A. R., Braun R., eds, 1999, Science with the Square 
Kilometer Array : a next generation world radio observa- 
tory 

Thompson A. R., Moran J. M., Swenson Jr. G. W., 2001, 

Interferometry and Synthesis in Radio Astronomy, 2nd 

Edition. Interferometry and synthesis in radio astronomy 

by A. Richard Thompson, James M. Moran, and George 

W. Swenson, Jr. 2nd ed. New York : Wiley, c2001.xxiii, 

692 p. : ill. ; 25 cm. "A Wiley-Interscience publication." 

Includes bibliographical references and indexes. ISBN : 

0471254924 
Verde et al. 2002, MNRAS, 335, 432 

Verheijen M. A. W., van Gorkom J. H., Szomoru A., 

Dwarakanath K. S., Poggianti B. M., Schiminovich D., 

2007, ArXiv e-prints, 708 
Zhang P., Pen U.-L., 2006, MNRAS, 367, 169 
Zwaan M. A., Meyer M. J., Staveley-Smith L., Webster 

R. L., 2005, MNRAS, 359, L30 
Zwaan M. A., van Dokkum P. G., Verheijen M. A. W., 

2001, Science, 293, 1800 
Zwaan et al. 2003, AJ, 125, 2842 



© 0000 RAS, MNRAS 000, 000-000 



